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In this paper, an M-analytical function and the sequence of its Taylor polynomials 
(which are Lyapunov functions different from those of Vanelli & Vidyasagar (1985, 
Automatica, 21(1):6 9-80)) is presented, in order to determine and approximate 
the domain of attraction of the exponentially asymptotically stable zero steady 
state of an autonomous, R-analytical system of differential equations. The analytical 
function and the sequence of its Taylor polynomials are constructed by recurrence 
formulae using the coefficients of the power series expansion of / at 0. 
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1 Introduction 



Let be the following system of differential equations: 
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where / : M" ^ R" is a function of class on W with /(O) = (i.e. a; = 
is a steady state of (1)). If the steady state x = is asymptotically stable [1], 
then the set Da{0) of all initial states for which the solution x{t;0,x^) of 
the initial value problem: 

X = f{x) x{0) = a;° (2) 



tends to as t tends to cxd, is open and connected and it is called the domain 
of attraction (domain of asymptotic stability [1]) of 0. 

The results of Barbashin [2], Barbashin-Krasovskii [3] and of Zubov ([4], The- 
orem 19, pp. 52-53, [5]), have probably been the first results concerning the 
exact determination of Da{0). In our context, the theorem of Zubov is the 
following: 

Theorem 1 An invariant and open set S containing the origin and included 

in the hypersphere B{r) = {x E M" : < r}, r > 0, coincides with the 
domain of attraction Da{0) if and only if there exist two functions V and ip 
with the following properties: 

(1) the function V is defined and continuous on S, and the function is 

defined and continuous on 

(2) -1 < V{x) < for any x e S\{0} and ip{x) > 0, for any x e M" \ {0} 

(3) \imV(x) = and Mmipix) = 

x^O x^O 

(4) for any 72 > small enough, there exist 71 > and cti > such that 
V{x) < —71 and ip{x) > ai, for \\x\\ > 72 

(5) for any y e dS, lim V{x) — —1 

(6) iV{x{t; 0; x')) = /(x(i; 0; x°))[l + V{x{t; 0; x'))] 

Remark 2 At this level of generality, the effective determination of Da{0) 
using the functions V and ip from Zubov 's theorem is not possible, because the 
function V (if ip is chosen) is constructed by the method of characteristics, 
using the solutions of system (1). This fact implicitly requests the knowledge 
of the domain of attraction Da(0) itself. 

Another interesting result concerning the exact determination of i?a(0), un- 
der the hypothesis that the real parts of the eigenvalues of the matrix |f (0) 
are negative, is due to Knobloch and Kappel [6]. In our context, Knobloch- 
Kappel's theorem is the following: 

Theorem 3 If the real parts of the eigenvalues of the matrix ||(0) are nega- 
tive, then for any function ^ : R" — > with the following properties: 

(1) (is of class onR"" 

(2) C(0) = and C,{x) > 0, for any x ^0 
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(3) the function ( has a positive lower limit on every subset of the set {x : 
\\x\\ >£},£> 

there exists a unique function V of class on -Da(O) which satisfies 

a. {VV{x),f(x))^-ax) 

b. V{0) = 

In addition, V satisfies the following conditions: 

c. V{x) > 0, for any x 7^ 

d. lim V{x) — 00, for any y e dDa{0) or for \\x\\ — > 00 

Remark 4 The effective determination of Da{0) using the functions V and 
( from Knobloch-Kappel's theorem (at this level of generality) is not possible, 
because the function V (if ( is chosen) is constructed by the method of char- 
acteristics using the solutions of system (1). This fact implicitly requests the 
knowledge of Da{0). 

Vanelli and Vidyasagar have established in [7] a result concerning the exis- 
tence of a maximal Lyapunov function (which characterizes Da{0)), and of 
a sequence of Lyapunov functions which can be used for approximating the 
domain of attraction Da{0). In the context of our paper, the theorem of Vanelli- 
Vidyasagar is the following: 

Theorem 5 An open set S which contains the origin coincides with the do- 
main of asymptotic stability of the asymptotically stable steady state x = 0, 
if and only if there exists a continuous function V : S ^ and a positive 
definite function ip on S with the following properties: 

(1) V{f)) = and V{x) > 0, for any x E S \ {0} (V is positive definite on S) 

(2) DrV{x'') = lim = _^(a;0)^ for any x" e S 

(3) lim V^(a;) = oo, for any y e dS or for \\x\\ 00 

Remark 6 The determination of Da{0) using the functions V and ip from 
Vanelli- Vidyasagar's theorem is not possible, for the same reason as in the 
case of the theorems of Zubov and Knobloch-Kappel. 

Remark 7 Restraining generality, and considering the case of an M.- analytic 
function f , for which the real parts of the eigenvalues of the matrix |f (0) are 
negative, Vanelli and Vidyasagar [7] establish a second theorem which provides 
a sequence of Lyapunov functions, which are not necessarily maximal, but can 
be used in order to approximate Da{0). These Lyapunov functions are of the 
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form: 

Tr. ^t"! -I- To ^t"! -I- -I- r.,, ( t\ 

(3) 



^ l + ?i(x) + g2(a:) + ... + g„»(a;) 



where ri and qi are i-th degree homogeneous polynomials, constructed using 
the elements of the matrix §^(0), of a positively definite matrix G and the 
nonlinear terms from the development of f . The algorithm of the construction 
of Vm is relatively complex, but does not suppose knowledge of the solutions of 
system (1). 

Very interesting results concerning the exact determination of the domains 
of attraction (asymptotic stabiUty domains) have been found by Gruyitch 
between 1985-1995. These resuhs can be found in [1], chap. 5. In these resuhs, 
the function V which characterizes the domain of attraction is constructed by 
the method of characteristics, which uses the solutions of system (1). Some 
illustrative examples are exceptions because for them V is found in a finite 
form, for some concrete functions /, but without a precise generally applicable 
rule. 

In the same year as VaneUi and Vidyasagar (1985), Bahnt [8], proved the 
following theorem: 

Theorem 8 ( see [8] or [9]) If the function f is R-analytic and the real parts of 
the eigenvalues of the matrix |^(0) are negative, then the domain of attraction 
Da{0) of the asymptotically stable steady state x — coincides with the natural 
domain of analyticity of the "M.- analytical function V defined by 

{VV{x)J{x)) = -\\x\\^ nO) = (4) 



The function V is strictly positive on Da{0) \ {0} and lim V{x) = oo for any 
y e dDa{0) or for \\x\\ —>■ oo. 

Remcirk 9 In the case when the matrix |f (0) is diagonalizable, recurrence 
formulae have been established in [10] (see also [9]) for the computation of the 
coefficients of the power series expansion in of the function V defined by (4) 
(called optimal Lyapunov function in [10]): 

Consider 5" : C" — > C" an isomorphism which reduces |f (0) to the diagonal 
form S-^%{0)S = diag{\i, A2...A„). Let be g ^ S'^ o f oS andW ^VoS. 
If the expansion ofW at is 

oo 

W{zi,Z2,...,Zn)=Yl Bnh-jA'^'-zt (5) 

m=2 \j\=m 
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and the expansions at of the scalar components Qi of g are 



gi{zi, Z2, Zn) — XiZi + X] X] ^}ij2.-jn^l^^2 

m=2 \j\=m 



(6) 



then the coefficients Bj-^j^ j^ of the development (5) are given by the following 
relations: 



i=l 



B 



3132— Jn 



E SipSiq if I j| = 2 and jp = = 1 



i=l 



lil-1 



Bji-kl...ji-ki+l...jn-kn] if \j\ > 3 



(7) 



Using these recurrence formulae, the optimal Lyapunov functions V and the 
domains of attraction Da{0) for some two-dimensional systems have been 
found in [10] in a finite form. 



Example 10 

xi = -Xxi + pix\ + P2X1X2 

X2 = —\X2 + P1X1X2 + P2X2 



A > 0,pi,pi e 



(8) 



The Lyapunov function corresponding to the zero asymptotically stable steady 
state of this system is 



V{xi,X2) 



rp^ I ^2 



In 



A 



A 



A {pixi + ^22:^2)^ A - {pixi + ^23:^2) Pixi + ^23:^2 



](9) 



and the domain of attraction is 

^a(O) = {X e : piXi + ^22^2 < A} 



(10) 



Example 11 

Xi — —Xxi + px\ + PX1X2 

X2 = —\X2 + px'{x2 + PX% 



A > 0,p e 
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The Lyapunov function corresponding to the zero asymptotically stable steady 
state of this system is 

V{xi, X2) = 77- In T 7^- — 2T (12) 

2p X — p[Xi +X2) 



and the domain of attraction is 

Da{0)^{xeR^:X-p{xl + xl)>0} (13) 



Therefore, when the function / is R-analytic, the real parts of the eigenvalues 
of the matrix |^(0) are negative, and the matrix |^(0) is diagonalizable, then 
the optimal Lyapunov function V can be found theoretically by computing the 
coefficients of its power series expansion at 0, without knowing the solutions 
of system (1). More precisely, in this way, the "embryo" Vq (i.e. the sum of the 
series) of the function V is found theoretically on the domain of convergence 
Dq of the power series expansion. A formula for determining the region of 
convergence Dq C -0^(0) of the series of V can be found in [11] or [9]. If Dq 
is a strict part of -Da(O), then the "embryo" Vq can be prolonged using the 
algorithm of prolongation of analytic functions: 

If Dq is strictly contained in Da{0), then there exists a point x^ G ODq such 
that the function Vq is bounded on a neighborhood of Let be a point x^ G 
Dq close to x^, and the power series development of Vq in Xi (the coefficients 
of this development are determined by the derivatives of Vq in x?). Using the 
formula from [11] or [9], the domain of convergence Di of the series centered 
in Xi is obtained, which gives a new part Di \ {DqCiDi) of the domain of 
attraction Da{0). The sum Vi of the series centered in x^ is a prolongation 
of the function Vq to Di and coincides with V on Di. At this step, the part 
Dq U Di of Da{0) and the restriction of F to £)o U Di are obtained. 

If there exists a point x^ G d{DQ\J Di) such that the function ^|d(,|J£)^ is 
bounded on a neighborhood of x^, then the domain Dq [j Di is strictly included 
in the domain of attraction i^a(O). In this case, the procedure described above 
is repeated, in a point x\ close to x^. 

The procedure cannot be continued in the case when it is found that on 
the boundary of the domain U -Di U ••• U -Dp obtained at step p, there are 
no points having neighborhoods on which V^|DoUDiy...yDp is bounded. We 
illustrate this process in the following example: 

Example 12 Consider the following differential equation: 

x^x{x-l){x + 2) (14) 
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a; = is an asymptotically stable steady state for this equation. The coeffi- 
cients of the power series development in of the optimal Lyapunov function 
are computed using (7): An = ^ sn^"-^^ ■ ^ 2. The domain of convergence 
Dq = (—1,1) of the series is found using the formula: 

xeDo iff ImI^|A„x"| < 1 (15) 

The embryo Vo{x) is unbounded in 1 and bounded in —1, as Vo(— 1) = 
We expand Vo{x) in —0.9 close to —1. The coefficients of the series centered in 
—0.9 are: A'„ = -^^[j^Yy^ + ^^r]- The domain of convergence Di of the series 
centered in —0.9 is given by: 

xeDi iff I5^^|A;(a; + 0.9)"| < 1 (16) 

and it is Di = (—2,0.2). So far, we have obtained the part D = DqIJDi — 
(—2, 1) of the domain of attraction -Da(O). As the function V is unbounded at 
both ends of the interval, we conclude that -Da(O) = (—2, 1). 

We have illustrated how this approximation technique described in [8,10,11] 
works in some particular cases. In more complex cases (for example if the 
right hand side terms in (1) are just polynomials of second degree), we can 
only compute effectively the coefficients Aj^^j^ j^ of the expansion of V up to 
a finite degree p. With these coefficients, the Taylor polynomial of degree p 
corresponding to V: 

p 

VQ{Xi,X2,...,Xn)= ^ ^ Aj^j^^^j^Xi^X2 ...X^j^ (17) 

m=2 \j\=m 

can be constructed. In the foUowings, it will be shown how Vq can be used in 
order to approximate Da{0). 



2 Theoretical results 

For r > 0, we denote by B{r) — {x E M.^ : \\x\\ < r} the hypersphere of radius 
r. 

Theorem 13 For any p > 2, there exists rp > such that for any x e 
B{rp) \ {0} one has: 

(1) Vp{x) > 

(2) {VVp{x),f{x)) <0 
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PROOF. First, we will prove that for p = 2, the function V2 has the proper- 
ties 1. and 2. For this, write the function / as: 

df 

f{x) = Ax + g{x) with A = -^(0) (18) 
and the equation 

{VV(x),f(x))^-\\x\\' (19) 

as 

{VV2{x),Ax) + (V(y - V2){x),Ax + g{x)) + {VV2{x), g{x)) = -\\x\\\20) 

Equating the terms of second degree, we obtain: 

(V^(x),Ax) = (21) 

As 1^2(0) = 0, it results that: 

00 

V2{x) = J \\e^'x\\^dt (22) 


This shows that ^(a;) > for any a; e \ {0}. 
On the other hand, one has: 

{WV2{x),f{x))^{VV2{x),Ax) + {VV2{x),g{x)) = 
^ -WxW^ + {VV2{x), g{x)) = 

^_||^||2[i_(VlM:^M^] (23) 

As hm i^v^i^lMx)) ^ Q ^Yiere exists > such that for any x e S(^\ {0}, 

11x11^0 ii^ii 

we have < \. Therefore, for any x e B{r2) \ {0}, wc get that: 



{^V2{x\m)<-h^xf (24) 

We will show that for any p > 2, the function satisfies conditions 1. and 2. 
Write the function as 

vM-y^{^)\^^^-^^:;^\ ^7^0 (25) 
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As lim ^''^^V2{x)^'^^ ~ ^' there exists such that for any x e B{r^) \ {0}, we 

Vp(x)-V2ix)i ^ 1 



have |li(£2_^| < i. Therefore, for any x e B{rl) \ {0}, we have 



V2{x) 

Vp{x) > W2{x) > (26) 



thus, Vp satisfies condition 1. 
On the other hand, we have: 



{WVp{x),fix))^{WV2{x),Ax)[l + 



{V{Vp-V2){x),f{x)) + {VV2{x),g{x))^ 



{VV2{x),Ax) 
{V{Vp - V2){x),f{x)) + {VV2{x),g{x))^ 



IFI, 

As hm = 0, there exists such that for any 

X e {0}, we have | (v(i^-i^)(.)./(x))+(vv^(.),g(.)) | ^ |_ Therefore, for any 

X e -B(r^) \ {0}, we have: 

(Vyp(x),/(x))<-^||x||^ (28) 



Therefore, for any x G B{rp) \ {0}, where Vp = m.m{rp,rp}, the function Vp 
satisfies conditions 1. and 2. 



Corollary 14 For an?/ p > 2, there exists a maximal domain Gp C suc/i 
that & Gp and for any x & Gp \ {0}, function Vp verifies 1. and 2. from 
Theorem 13. In other words, for any p > 2 the function Vp is a Lyapunov 
function for (1) (in the sense of [1]). 

Remcirk 15 Theorem 13 provides that the Taylor polynomials of degree p > 
2 associated to V in are Lyapunov functions. This sequence of Lyapunov 
functions is different of that provided by Vanelli and Vidyasagar in [7]. 

Theorem 16 For any p > 2, there exists c > and a closed and connected 
set S of points from x e IR", with the following properties: 

(1) Oe Int{S) 

(2) Vp{x) < c for any x G Int{S) 

(3) Vp{x) = c for any x E dS 

(4) S is compact and included in the set Gp. 
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PROOF. Let be p > 2 and Vp > determined in Theorem 13. Let be c = 
min Vp{x) and S' — {x E B{rp) : Vp{x) < c}. It is obvious that c > and 

\\x\\=rp 

that there exist x* with = such that V{x*) — c. The set S' is open, 
OeS' and S' C 'B{^ C Gp. 

We will prove that Vp{x) — c for any x e dS'. Let be x e 95". Thus, ||x|| < Vp 
and there exists a sequence x'^ G S' such that a;*^ — >• a; as /c — > oo. As lp(a;'^) < c, 
we have that Vp{x) — lim V^(a;*^) < c. The case = Vp and Vp(a;) < c is 

fe— iOO 

impossible, because c = min Vp{x). The case ||x|| < and V^(x) < c is also 

ll^ll=''p 

impossible, because this would mean that x belongs to the interior of the set 
S', and not to its boundary. Therefore, for any x e dS' we have Vp{x) — c. 

If the set S' is not connected (see Example 29 in this paper), we denote by 5"' 
its connected component which contains the origin, and let he S — S". Then it 
is obvious that S is connected (being the closure of the open connected set S"), 
G Int{S) = S", and that for any x G Int{S) = S", we have Vp{x) < c. More, 
as dS = dS", we have Vp{x) = c for any x G dS. As S" is bounded, we obtain 
that the closed set S is also bounded, thus, it is compact. As S" C B{rp) C Gp, 
we have that S — S" C B{rp) C Gp. Therefore, S satisfies the properties 1-4. 

Lemma 17 Let be p > 2, c > and a closed and connected set S satisfying 
1-4 from Theorem 16. Then for any x^ G S, the solution x{t; 0,x^) of system 
(1) starting from x^ is defined on [0, oo) and belongs to Int{S) for any t > 0. 



PROOF. Let be x^ G S. We denote by [O,p,j.o) the right maximal interval of 
existence of the solution x{t; 0, x^) of system (1) with starting state 

First, if x° G /nt(5)\{0}, we show that a;(t;0,a;°) G Int{S), for alH G [Q^P^o)- 
Suppose the contrary, i.e. there exists T G (0,/3j.o) such that x{t;0,x^) G 
Int{S), for t e [0,T) and x(T;0,x°) G dS (i.e. Vp{x{T;0, x^)) = c). As 
x(t;0,x'^) G Gp \ {0}, for t G [0,T), Vp{x(t;0,x^)) is strictly decreasing, 
and it follows that Vp{x{t;Q, x^)) < Vp{x^) < c, for t G (0,T). Therefore 
Vp{x{T;0,x'^)) < c, which contradicts the supposition a;(T;0,a;°) G dS. Thus, 
x{t; 0, x^) G Int{S), for all i G [0, /3^o). (It is clear that for x° = 0, the solution 
x{t; 0, 0) = G Int{S), for all t > 0.) 

If x^ G 95, we show that x{t;0,x'^) G /nt(5'), for all t G (0,/5,,o). As the 
compact set 5 is a subset of the domain Gp, the continuity of x{t; 0, 
provides the existence of Tj.o > such that x{t; 0, G Gp \ {0} for any 
t G [0,T^()] c [0,j3xo). Therefore Vp{x{t;0,x^)) is strictly decreasing on [0,Tj.o], 
and it follows that Vp{x{t;0,x^)) < lp(x°) = c, for any t G (0,Tj,()). This 
means that Vp{x{t;0,x^)) G Int{S), for any t G (0,Ta.o]. The first part of the 
proof guarantees that x{t;0,x^) G Int{S), for all t G [T!j,o, /3^o), therefore, for 
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all t e (0,/3^o). 



In conclusion, for any G S, we have that x{t;0,x'^) G Int{S), for all t G 
(0,/3.o). 

As for any x^ G S*, the solution x{t;0,x^) defined on [0,/9a.o), belongs to the 
compact S, we obtain that /S^-o = oo and the solution x{t]0,x'^) is defined on 
[0, oo), for each x^ G S. More, x{t; 0, x°) G /nt(5), for all t > 0. 

Remcirk 18 Lemma 17 states that a closed and connected set S satisfying 1-4 
from Theorem 16 is positively invariant to the flow of system (1). 

Theorem 19 (LaSalle-type theorem) Let be p > 2, c > and a closed and 
connected set S satisfying 1-4 from Theorem 16. Then S is a part of the of 
the domain of attraction Da(0). 

PROOF. Let be x° G 5'\{0}. To prove that lim x{t\ 0, x^) = 0, it is sufficient 

t — ^oo 

to prove that lim x{tk;0,x^) = 0, for any sequence — > oo. 

k—*oo 

Consider tk — > oo. The terms of the sequence x{tk; 0, belong to the compact 
S. Thus, there exits a convergent subsequence x{tkj ; 0, x°) — > y° G -S. 

It can be shown that 

Vp{x(t; 0, > Vp(y'^) for alH > (29) 

For this, observe that x{tk/,Q, x^) y^ and Vp is strictly decreasing along the 
trajectories, which implies that Vp{x{tkj'iO,x'^)) > Vp{y^) for any kj. On the 
other hand, for any t > 0, there exists kj such that tk > t, and therefore 
Vp{x{t;0,x^)) > Vp{x{tk,;0,x^)) > i;(/). 

We show now that y° = 0. Suppose the contrary, i.e. y° 7^ 0. Inequality (29) 
becomes 

Vp{x{t; 0, a;°)) > Vp{y°) > for aU t > (30) 

As Vp{x{s; 0, y^)) is strictly decreasing on [0, 00), we find that 

Vp{x{s; 0, 1/°)) < V;(/) for all s > (31) 

For s > 0, there exists a neighborhood ?7x(s;0,?;0) C 5 of a;(s;0,|/°) such that 
for any x G Ux(s]0,yO) we have < Vp{x) < Vp{y'^). On the other hand, for the 
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neighborhood Ux{s;0.y") there exists a neighborhood UyO G S oi such that 
x{s; 0, y) e U^{s;o,yO) for any y eUyO. Therefore: 

Vp{x{s; 0, y)) < Vp{y'') for all y E UyO (32) 

As x(tkj;0,x^) — > |/°, there exists kj such that a;(tfc^.; 0, a;°) e t/j^o, for any 
/cj > kj. Making y — ^(tfc^.; 0, x°) in (32), it results that 

Vp{x{s + tk-.;0, x')) = Vp{x{s; 0, x{tk.; 0, x°))) < i;(y°) for kj > kj{33) 

which contradicts (30). This means that = 0, consequently, every conver- 
gent subsequence of a;(tfe; 0, a;°) converges to 0. This provides that the sequence 
x(tk;0,x^) is convergent to 0, for any — >■ oo, thus lim x{t;0,x^) = 0, and 

t— »oo 

Therefore, the set S is contained in the domain of attraction of Da{0). 

Corollary 20 For any p > 2 and c > there exists at most one closed and 
connected set satisfying 1-4 from Theorem 16. 

PROOF. Suppose the contrary, i.e. for a p > 2 and c > there exist two 
different closed and connected sets 5*1 and 5*2 satisfying 1-4 from Theorem 
16. Assume for example that there exists x^ E Si \ 82- Due to Theorem 19, 
Si C -Da(O) and therefore lim x{t] 0, a;°) = 0. As x^ ^ 5*2, and 5*2 is a closed and 

t— »oo 

connected neighborhood of 0, there exists T > such that a;(T;0,a;*') G 882. 
Therefore, Vp{x{T;{),x^)) — c which contradicts Lemma 17. Consequently, we 
have Si C 82- By the same reasons, 82 Q 81. Finally, 81 — 82. 

Remark 21 If for p > 2 and c > there exists a closed and connected set 
satisfying 1-4 from Theorem 16, then it is unique and it will be denoted by N^. 

Corollary 22 Any set is included in the domain of attraction Da{0). 

Lemma 23 Let be p > 2 and c > such that there exists the set N^. Then, 
for any d e (0, c] the set {x G iV^ : Vp{x) < c'} coincides with the set . 

PROOF. Let be d G (0,c]. It is obvious that N'^ is included in the set 
{x G iVp" : Vp{x) < d}. Let be x^ G such that Vp{x^) < d. We know that 
Vp{x{t;0,x^)) < Vp{x'^) < d, for any t > 0. Theorem 19 provides that x^ G 
Np C -Da(O), therefore, x° is connected to through the continuous trajectory 
x{t; 0, x°), along which Vp takes values below d. In conclusion, x° G N^' . 
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Theorem 24 If for p > 2 and c > there exists N^, then for any c' e (0, c) 
there exists N^' and N^' C N^. More, for any Ci,C2 G (0, c) we have N^^ C N^'^ 
if and only if ci < c^- 



PROOF. Lemma 23 provides that for any d e (0, c) there exists — {x ^ 
■■ Vp{x) < c'}. It is obvious that C N^. 

Let's show that for any ci, C2 G (0, c) we have N^^ C N^^ if and only if ci < C2. 

To show the necessity, let's suppose the contrary, i.e. N'^'^ C A'^^ and ci > C2. 
Let be x^ e SA/"^^ c Gp. Then = C2 and as x° e Gp, we get that 

Vp{x{t] 0, < V;,(a;°) = C2 for any t > (34) 



Theorem 19 provides that x^ G dN'^^ C Da(0), therefore lim a;(t;0,x°) = 0. 
As A'pi and A'^^ are connected neighborhoods of and A'^^ C N^^ , there exists 
T > such that x{T- 0, x°) e aA"^i . This means that Vp{x{T- 0, a;°)) = ci > C2, 
and (34) provides that ci = C2. As A^^ is strictly included in A^^^ there exists 
X G c^A'pi (i.e. Vp{x) = c\ = C'l) such that x G Int{Np^). This contradicts the 
property 2 from Theorem 16 concerning A^^ conclusion, Ci < C2- 

To prove the sufficiency, let's suppose that ci < C2 and let be x^ G A^^ \ {0}. 
As x° G Afi C Da(0), we have that hm x(t; 0, x°) = 0, so x° is connected to 

^ t— >oo 

through the continuous trajectory x{t;0,x^). More, as x^ G A^^ \ {0}, we 
have Vplx'^) < Ci < c-^. This means that ,x° G A'^^, therefore A^^ C N'^ . The 
inclusion is strict, because A'^^ = A'^^ means dN^ = dN^^ , i.e. Ci = C2, which 
contradicts ci < C2. 



Corollary 25 For a ^'wen p > 2, the set of all N^-s is totally ordered and 
[JNp is included in Da,{0). Therefore, for a given p > 2, the largest part of 

c 

Da{0) which can be found by this method is U A^- 

c 



For any p > 2 let be i?p = {r > : B{r) C Gp}. For r G -Rp we denote by 



Corollary 26 For any r G Rp, there exists the set Np^ and Np^ C B{r). 

Corollary 27 For any p > 2 and any r' , r" G Rp, r' < r" , such that Vp is 
radially increasing on B{r") we have < . 

Remark 28 In some cases, it can he shown that the function Vp is radially 
increasing on Gp.- 
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a. V2 is radially increasing on M."'; 

b. If n — 1, then for any p > 2, Vp is radially increasing on Gp. 
This result is not true in general, provided by the following example: 
Example 29 Let be the following system of differential equations: 



Xi = -Xi - X1X2 

(35) 

X2 = -X2 + X1X2 



for which (0, 0) is an asymptotically stable steady state. For p — 3 the Lya- 
punov function \3{xi,X2) is given by: 

1 1 

V3{xi,X2) = 2 (^1 + ^2) + 3(^1^2 - ^2x1) (36) 



Consider the point (3-\/5, e dGs and let be : [0, 1) — > Gs defined by 
g{X) = VsiSVEX, V5X). The function g is increasing on [0, ^] and decreasing 
on (^,1), therefore, the Lyapunov function V3 is not radially increasing on 
the direction (3-\/5, -\/5)- In conclusion, V3 is not radially increasing on 6*3. 

More, for c = 0.32 there exists N^, but the set {x = {xi, X2) G G3 : V3{xi, X2) < 
c} is not connected. The reason is that the point {xi,X2) — (^, H) G dGs 
with V3{xi,X2) = has a nonempty neighborhood U such that V3{xi,X2) < c, 
for any (xi, X2) e G^^U and (Gg nU)n = 0. 

Theorem 30 For any p > 2 there exists pp > such that Vp is radially 
increasing on B{pp). 



PROOF. It can be easily verified that V2 is radially increasing on R", using 
relation (22). This provides that for any x G M" \ {0}, the function g2 : 
R+ R+ defined by fiff (A) = V2{Xx) is strictly increasing on R+, therefore 
^^f (A) > on R* , i.e. 

{W2{Xx),x) > for any A > and a; G R" \ {0} (37) 



Let be p > 2, x G R" \ {0} and g^:R+^ R+ defined by g^{X) = Vp{Xx). One 
has: 



^9;{X) = {VVp{Xx),x) = {VV2{Xx),x) + {V{Vp - V2){Xx),x) 
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= {VV,{Xx),x)il + ^^y^^^^^^^^ ) (38) 

As lim = 0, there exists p, > such that < |, 

for any x e -B(pp) \ {0}. Relation (38) provides that for any x e B{pp) \ {0}, 
we have: 

^^^(A) > ^{WV2{Xx),x) > for any A > (39) 



Therefore, for any x e B{pp) \ {0}, the function is strictly increasing on 
R", i.e. Vp is radially increasing on B{pp). 

Theorem 31 Let be p > 2 and c > such that there exists the set N^. 
Suppose that for any d < c, the sets N^' have the star-property, i.e. for any 
X G Np and for any X G [0, 1) one has Xx G Int{Np'). Then Vp is radially 
increasing on N^. 



PROOF. Let be a;° G dN^ and < Ai < Ag < 1. We have to show that 
Vp{Xix^) < Vp{X2X°). Denote ci = Vp{Xix^) > 0, C2 = Vp{X2X°) > and 
suppose the contrary, i.e. ci > C2. Theorem 24 provides that N^^ D N^^. 
Lemma 23 guarantees that X2X G dN^^ . As N'^'^ has the star- property, then for 
A = ^ G (0, 1), we have that A(A2a:) = XiX G Int^N^^), so Ci = Vp{Xix) < C2 
which contradicts the supposition ci > C2. Therefore, Vp is radially increasing 
on N^. 

Remcirk 32 a. For any x G Dq, there exists > 2 such that x & Gp for any 

V > Px,- 

b. If n = 1, there exists Pq>2 such that Dq C Gp, for any p > po- 

c. If there exists r > such that B{r) C Gp for any p > 2, then there exists 
Po>2 such that Dq C Gp^. 

Conjecture 33 For any x G Da{0) there exists p > 2 and c > such that 



3 Numerical example: the Van der Pol system 



We consider the following system of differential equations: 

Xi = -X2 

X2 ^ Xi- X2+ x\x2 



(40) 
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The (0, 0) steady state of (40) is asymptotically stable. The boundary of the 
domain of attraction of (0, 0) is a limit cycle of (40). 

For p = 20 we have computed that the largest value c > for which there 
exists the set A^^^ is C20 = 8.8466. For p = 50, the largest value c > for which 
there exists the set is C50 = 13.887. In the figures below, the thick black 
curve represents the boundary of Da{0, 0), the thin black curve represents the 
boundary of Gp and the gray surface represents the set N^p. The set Np^° 
approximates very well the domain of attraction of (0, 0). 




-3-2-10123 -3-2-10123 



Fig. 1. The sets N^^", G20 and Da{0, 0) for Fig. 2. The sets N^^'', G50 and £)„(0, 0) for 
system (40) system (40) 
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